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Nomenclature 


Mathematical Constants 


c 7 

= 2e -7 

1.122919... 

e 

base of natural logarithms 

2.7182818... 

7 

Euler’s constant 

0.577215665... 

7 r 

circumference /diameter ratio 

3.14159265... 

Physical Constants 


c 

speed of light in vacuum 

2.99792458 x 10 10 cm/s 

e 

electric charge 

1.602176487 x 10" 19 C 

h 

Planck constant 

6.6261 x 10~ 2 ' erg s 

h 

Planck constant, reduced, =h/2n 

1.0546 x 10~ 2 ' erg s 

rn e 

electron rest mass 

9.1094 x 10" 28 g 

N 0 

Avogadro’s constant 

6.0221 x 10 23 particles/mol 

r e 

classical electron radius 

2.8179 x 10“ 13 cm 

a 

fine structure constant 

7.2974 x 10 -3 


General Variables Note: if units are not specified, then the quantity is unitless. 

A atomic mass (g/mole) 

Ap curve-fit parameter (MeV -p ) 

A pp pair production curve-fit parameter (cm 2 /g) 

j4g annihilation factor in positron production term 

Ap positron annihilation loss term 

A w curve-fit parameter used in energy straggling (MeV) 

a, b redefinition variables 

a' energy parameter 

a, moments of scattering power for electron transport (g/crn 2 )' 

ap parameter used to define the density effect correction to the 

Bethe-Bloch formula 

C term in modified Bethe-Bloch formula 

Ce parameter used to define applicability of density correction term to 

Bethe-Bloch formula 

D absorbed radiation dose (Sv) 

di depth in the material slab of the ith point (g/cm 2 ) 

E total energy (kinetic + rest mass) (MeV) 

E p material plasma energy (MeV) 

E 4 =4 MeV 

F cutoff dependent correction to Bethe-Bloch formula 

F c , f c , fee Coulomb correction terms to bremsstrahlung cross section 

/ fraction of total photon attenuation for a given process 

/i , / 2 , fs energy functions used in the definition of the cutoff correction F 
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5' 

e 

c 

v 

e 

O r 

K 

A 

Atr 


ratio of energy terms in Bethe-Bloch formula 

energy function for bremsstrahlung production 

energy parameters used in the annihilation cross section 

mean atomic ionization potential (MeV) 

integration constant 

curve-fit parameter 

triplet term parameters defined in Table 1 

slope parameter for photoelectric absorption coefficient 

number density (particles/cm 3 ) 

particles per unit mass (particles/g) 

electron mass scattering power (cm 2 /g) 

curve-fit parameter 

functions in pair production cross section formula (MeV 2 ) 

ratio of photon energy to electron rest mass energy 

range of electron in given material (g/cm 2 ) 

scaled energy ratio 

electron deflection radius (g/cm 2 ) 

stopping power (MeV cm 2 /g) 

curvilinear distance traversed by an electron in a material (g/cm 2 ) 

particle kinetic energy (MeV) 

particle speed (cm/s) 

residual energy (MeV) 

scaled thickness of a material slab (g/cm 2 ) 

position or distance (g/cm 2 ) 

scaled energy variable 

scattering function 

atomic number 

axial (projected) distance traversed (g/cm 2 ) 
ratio of particle speed to the speed of light in vacuum 
ratio of particle total energy to rest mass energy 
energy parameter (MeV) 

density correction term in Bethe-Bloch equation 

energy function used in determining screening functions used in bremsstrahlung 
cross section 

density correction term in pair production cross section 
screening function in Rutherford cross section 
source term (particles/ (MeV g)) 
transmission function 

scattering angle; also, polar angle in spherical coordinates 
scattering angle ratio 

ratio of photon energy to electron rest mass energy 
shell correction term in Bethe-Bloch formula 
transport mean free path (g/cm 2 ) 


IV 



p 

photon energy deposition coefficient (cm 2 /g) 

e 

bremsstrahlung triplet production variable (cm 2 /g) 

p 

material density (cm -3 ) 

p' 

scaled mass density (cm 6 /g 2 ) 

PS 

ratio of stopping powers 

o 

cross section (cm 2 /g) 

T 

ratio of kinetic energy to rest mass energy 

<f* 

Klein-Nishina correction function 


flux (particles/ (cm 2 MeV)) 

0 

spherical coordinate azimuth angle 

01) 02 

screening functions 

UJ 

electron orbital frequency (s -1 ) 

N 

deviation of transport mean free path (g/cm 2 ) 

Subscripts 


ann 

refers to a positron annihilation process 

B 

reference to Bohr stopping power formulation 

CSDA 

continuous slowing down approximation 

col 

indicates collision process 

e 

refers to electrons 

en 

refers to energy deposition in photon transport 

I 

reference to average ionization potential 

in 

refers to an incoherent scattering process 

IC 

refers to atomic K shell 

KN 

refers to Klein-Nishina formula 

L 

refers to atomic L-shell 

1 

reference to lesser (or lower) value 

max 

maximum allowed value 

min 

minimum allowed value 

PE 

refers to photoelectric process 

PP 

refers to pair production process 

R 

refers to Rutherford formula 

rad 

denotes radiative (bremsstrahlung) process 

Tot 

indicates summation of referenced quantities 

ti- 

reference to transport process 

ll 

reference to higher (upper) value 

V 

refers to photons 

+ 

reference to positron 

Superscripts 


(b) 

indicates source quantity 

V 



(t) 

+r,± 
/ // 


indicates transmitted quantity 

reference to positron, electron, or both (respectively) 
denotes dummy integration variables 
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Abstract 


A deterministic computational procedure has been developed to describe transport of 
space environment electrons in various shield media. This code is an upgrade and 
extension of an earlier electron code. Whereas the former code was formulated on 
the basis of parametric functions derived from limited laboratory data, the present 
code utilizes well established theoretical representations to describe the relevant in- 
teractions and transport processes. The shield material specification has been made 
more general, as have the pertinent cross sections. A combined mean free path and 
average trajectory approach has been used in the transport formalism. Comparisons 
with Monte Carlo calculations are presented. 

1 Introduction 

The concept of developing a rapid analysis electron transport code at NASA Langley Research 
Center (LaRC) arose from a desire to have a companion code for the LaRC deterministic ion 
code HZETRN [1], The initial effort resulted in a code [2] incorporating parameterizations of 
the effective range of electrons. These parametric functions were derived from experimental 
range data for electron beams incident on a variety of materials. In addition, interpolation and 
extrapolation methods were used in the transport computational process. Although numerous 
comparison calculations indicated that the code was reasonably accurate for low Earth orbit 
(LEO) conditions, the limit of applicability was dictated by the data range on which the code 
algorithms were based. Applications to gas giant planet missions with very high energy trapped 
electron environments required questionable extrapolation with regard to both energy and shield 
composition. Thus, work on a substantial upgrade began and has resulted in the new formulation 
described herein. 

In the new code, the parametric effective range approach has been replaced in favor of a 
continuous slowing down approximation (CSDA) combined with an elastic multiple scattering 
formulation to define an electron mean free path and a transmission function at a given target 
location. The cross sections relevant to CSDA and multiple scattering are described in the 
following section. The pertinent cross sections are calculated with well established theoretical 
models, valid over a much greater energy range than were the expressions in the earlier code. 
Consequently, a calculation of positron production and annihilation has been incorporated. 
The cross section equations are cast in terms of specific elemental atomic species, from which 
cross sections applicable to any user defined molecular system are automatically constructed. 
Maximum accuracy would be obtained with an attachment of tabular cross section data, which 
are available from NIST or other sources. Databases inclusive of all elements would be inherently 
large and require coding to enable appropriate table look-up procedures. For the present version, 
the sacrifice of some precision for the sake of speed is warranted. Future versions of the code 
may be equipped with options to include tabular cross section data as desired by the user. 

Details of the transport formulation are described in a subsequent section and are based on 
motion related to some initial direction (or axis of propagation). Behavior of the slowing and 
stopping of electrons and their associated bremsstrahlung is evaluated in terms of quantities 
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contained in the direction of the hemisphere centered on the initial direction. The scattering 
and mean free path calculations permit definition of mean trajectories, relative to the initial 
motion direction, so that axisymmetric spatial and energy distributions may be inferred. 

Comparisons with corresponding Monte Carlo calculations are shown. The Monte Carlo 
cases utilized several hours of machine time, whereas the LaRC calculations were practically 
instantaneous. Despite some difference in final magnitudes, general functional behavior was 
consistent in the results, which means the use of the deterministic code in trade studies could be 
performed quickly and credibly. A further improvement in the present code is that the required 
input has been greatly simplified. The present code includes an atomic database permanent 
hie and requires two user defined hies specifying material composition along with the incident 
environment differential energy spectra for electrons and photons. Most results are presented for 
the selected materials: aluminum (Al), tantalum (Ta), and water (H2O). This selection of target 
materials was chosen to provide examples of light and heavy materials in addition to elemental 
and compound materials. The materials were also chosen for their uses in space applications. 
Aluminum is a common spacecraft structural material, water closely approximates human tissue, 
and tantalum is used in shielding sensitive electronics from electron and photon exposure in the 
planetary trapped radiation belts. 


2 Cross Sections 

The theoretical description of the propagation of energetic electrons and photons in condensed 
media ordinarily requires some representation of interaction cross sections (or associated pa- 
rameters) with constituent atoms of the medium. Application of the principles of quantum and 
semi-classical physics have gone far toward providing formulations for electron and photon in- 
teraction cross sections. Although many of these formulations are mathematically complicated, 
they are often amenable to approximation and parameterizations that greatly simplify their 
practical application. The following subsections describe the formulations used for interactions 
of electrons and photons. 

2.1 Electrons 

Free energetic electrons passing through a material are slowed and may eventually be stopped 
as they interact with the electric fields of nuclei and the bound electrons of the material. In- 
teractions also result in a change of direction (scattering). In the present analysis, the electron 
deceleration process is expressed in terms of energy loss by collisions and accompanying photon 
production (bremsstrahlung). The dominant scattering process is assumed to be that of small 
angle elastic (Rutherford) scattering. It is the only scattering process taken into account in the 
present formulation. 

2.1.1 Collisional Slowing and Stopping 

Electron deceleration by collision processes, in which energy from the projectile electron is 
imparted to the electrons of the medium, is usually cast in terms of the stopping power, S, 
which is the energy loss per unit scaled distance, —dE/dx. Note that the distance x is measured 
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in units of g/cm 2 . The resultant expression has been often referred to as the modified Bethe- 
Bloch formula and is given here in the notation of Ref. [3] , 

jp 

S'coi = = 2^r 2 e m e c 2 p\N 0 Z/A)G T C ± , (1) 

where 


and 


(T + m e c 2 ) 2 
T(T + 2 m e c 2 ) 





2 (T + 2 m e c 2 ) 
m e c 2 (I / m e c 2 ) 2 


+ F ± 


5 - X> . 


(2) 


The superscript (±) indicates applicability to both electrons (-) and positrons (+). Here, I is 
the mean atomic ionization potential, m e is the rest mass of the electron, c is the speed of light 
in vacuum, p' is the scaled mass density with distance measured in g/cm 2 , Nq is Avogadro’s 
constant, Z is the atomic number of the material, A is the atomic mass of the material, r e is 
the classical electron radius, T is the kinetic energy of the electron or positron, 5 is the density 
correction term, and A is the shell correction term. The F± quantity in Eq. (2) is a function 
only of projectile energy and is slightly different (a few percent) for electrons and positrons. 
The present formulation implements only the F± function that applies to electrons. It may be 
written as the sum of three terms, 


f ± « F- = h + h - h. 


( 3 ) 


where 


fi 


h 


h 


— — h In A( ^ ~ A ) ] 

T - A + L (rn e c 2 ) 2 \ ’ 

A 2 /2 + m e c 2 (2T + m e c 2 ) ln(l — A /T) 

(: T + m e c 2 ) 2 

T(T + 2i n e c 2 ) 

( T + m e c 2 ) 2 


and 


A = T/2. 

The quantity 5 in Eq. (2) represents a modification to the Bethe-Bloch stopping power usually 
referred to as the “density correction” and is presented in detail in Ref. [4], The formulation of 
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this correction term is based on extensive experimental data and is highly parameterized. The 
energy conditional relationships are cast in terms of an energy variable, 


V = logioM 


ln(7 2 l 2 ) 

2 ln( 10 ) 


In [(T 2 + 2 Tm e c 2 )/(m e c 2 ) 2 
2 ln( 10 ) 


(4) 


where (3 is the ratio of the electron speed to the speed of light and 7 = 1 + T/m e c 2 . The 
elemental “plasma energy,” E p , also enters into the formulation. Its value is given by 


E p = hc^/ATrr e p(N 0 Z/A) (5) 

and appears in the parameter Ce • Here, h is the reduced Planck constant and p is the material 
density. The parameter Ce is given as 


Ce = 1 + 2 ln(// E p ), ( 6 ) 

where I is the mean atomic ionization potential. The Ce parameter is used in conjunction 
with other parameters to define limits of application to the 5 correction term. The value of 5 is 
defined by 


6 = 


0 y <yo 

4.605l7y - C E + a H {yi - y) 3 yo < y < yi 
4.60517y - C E yi < y 


(7) 


and applies for electron energies > 0.0001 MeV (or 100 eV). The parameters yo, yi, and an are 
defined as 


= f 0.2 C E < 5.215 

yo = \ 0.326 C E - 1.5 C E > 5.215 ’ 

yi = 3.0, 

a H = 4.60517(y a - yo)/{yi - yo) 3 , 
y a = C E / 21 n( 10 ) = C s /4.60517. 


The remaining correction term, A, in Eq. (2), is called the “shell correction” and is only 
important at low energies when the projectile electron speed, v, approaches the speed of a 
bound electron, vj. Its value, as used in the present calculations, is given by [5] 


, ( v i\ 2 ^ 7(7 + 2 m e c 2 ){T + m e c 2 ) 2 

~U ~ T{T + 2m e c 2 ){I + m e c 2 ) 2 ' U 

It is seen in Eq. (1) that 5 co i is unbounded as T — > 0 due to the Gt factor. An early 
non-relativistic stopping power derivation was provided by Bohr [5], which approached zero 
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at very low energies and exhibited a maximum. The Bohr formula may be written in CGS 
(centimeters-grams-seconds) units as 


'S'Bohr 


47rne 4 

m e v 2 


In 




with Cj = 2e 7 and ui = I /h. From d5 j^ hr = 0, 


From Eq. (9), 


( 9 ) 


b = 2 3 / 2 Cry/(e 2 LUiy/rn^) 

and 


_ /e 3 / 2 eVm^A 2/3 

B \ 2‘ i / 2 Cryh ) ’ 

where T B is the projectile kinetic energy at peak stopping power. If Eq. (1) is compared with 
Eq. (9), the natural logarithm term in Eq. (9) can be identified with the C± in Eq. (1). This 
gives 


C± ax = In (&Tjf) , 

Sjnax = 2-Krlm e c 2 N 0 G T {Z/A)C^ x . 

For T < T B , an extension to low energies is provided by Sigmund [5]. A simple functional fit 
has been derived for the present work, 

Scoi = S max exp {-0.179 [In ( T B /T )] 2m } for T <T B (10) 

and Eq. (1) is applied for T >T B . 

Three simple materials important to space exploration applications have been selected to 
demonstrate results of the present model calculations: water, aluminum, and tantalum. They 
were chosen to represent a broad range of atomic properties and material densities. Fig. 1 depicts 
the collision stopping powers as calculated for these substances over the energy range of the code 
formulation. The most notable features are the generally higher stopping powers for lighter 
materials and the shift of peak stopping power to higher energy for the higher atomic number 
materials. For comparison, the stopping powers obtained from the database tabulations of the 
National Institute of Standards and Technology (NIST) [6] are also shown. The comparisons 
indicate that the present formulation is quite acceptable over the energy range of interest for 
space applications (0.001 - 1000 MeV), with the greatest discrepancies occurring at the lowest 
energies for heavier elements. 
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Electron Kinetic Energy, MeV 


Figure 1: Collisional stopping power from current model and NIST calculations [6]. 


2.1.2 Radiative Energy Loss (Bremsstrahlung) 

Accelerating (and decelerating) charged particles lose energy by photon emission. Elaborate 
quantum mechanical calculations have gone far toward quantifying bremsstrahlung cross sec- 
tions. A detailed description of the theoretical results is given in Ref. [7]. The formulation may 
be greatly simplified by using several parameterizations; the implementation for this work is 
essentially the same as that for the default cross sections of the EGSnrc Monte Carlo code [8]. 
Several pertinent parameters that lead to the ultimate cross section formula may be defined as 


where 


F, = 4 


3 HZ) + fc 


( 11 ) 


with 


f 0.0 E u < 50 MeV 
1 fee K > 50 MeV ’ 


f cc = (aZ) 2 J2 


n{n 2 + a 2 Z 2 ) 


“ 7.4564 x 10- 5 Z L9137 . 
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Here, a is the fine structure constant. The term F c serves to provide a high energy Coulomb 
correction. Another parameter defines the effects of screening in terms of an energy variable, 
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S 


/!36\ 

E u m e c 2 

UW 

_E e (E e — E u ) _ 


which involves the electron total energy, E e = T + m e c 2 , and the emitted photon energy, E v . 
The screening functions have been specified as [8] 


<f>i(S s ) 

4>2 (5 a ) 


20.867 - 3.2425s + 0.625#? 5 S < 1 
21.12- 4.184 ln(# + 0.952) S 9 > 1 

20.209 - 1.9355s + 0.086#? 5 S < 1 
<j>i(S s ) S s > 1 


( 12 ) 


A final parameter, # which accounts for the triplet production process, depends upon atomic 
charge and involves the Coulomb correction above, along with the quantities L and L' that have 
been tabulated in Ref. [8] and are given in Table 1 . A parameterized fit for £ has been developed 
for this work, 


i = L'/{L~f c ) 


1.147 sin(0.42 In Z) + 0.12 j [sin(0.76 In 



(13) 


The final expression for the cross section for production of a photon of energy E u by inter- 
action of an electron of total energy E e on an atom of charge Z may be written as 


where 


d(Ju,e _ «#? N 0 Z(Z + 0 h 
dE u E v A 


(14) 


H v = 


1 + 


E e — E u 

Ee 




(#2 - F c 


An obvious energy constraint is that only photons can be produced that have energy less 
than the initial electron kinetic energy. Consequently, da Uye /dE v is calculated at each electron 
energy T, for selected values of E u /T\ [0.001, 0.05, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, 0.9, 0.95, 0.98, 
0.99, 1.0]. This array has been chosen because the cross section variations are greater at either 
end of the energy range. A source term (, for production of photons by an electron flux ip e , may 
be written in terms of the bremsstrahlung cross section as 


T 

C (Eu,T) = J VeiT^^dT' for E v < T. (15) 

E u 
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Table 1: Parameter values used in Eq. (13) from Ref. [8]. 


z 

L 

L' 

1 

5.31 

6.144 

2 

4.79 

5.621 

3 

4.74 

5.805 

4 

4.71 

5.924 

> 4 

ln(184.15Z~ 1 /3) 

ln(1194Z~ 2 / 3 ) 


The corresponding energy loss term due to radiative processes may be evaluated as 


‘S'rad — 


' E'^dE' 


dE' v 


(16) 


Radiative loss stopping powers for the three materials selected for illustration are shown in 
Fig. 2. In contrast to the collision stopping powers, bremsstrahlung stopping power is greater 
for higher charge elements and increases monotonically with energy. The logarithmic stopping 
power for the present approximate formulation is very nearly linear with the logarithm of the 
energy and begins to have substantial effects only for kinetic energies greater than a few MeV. 
Comparison with NIST data reveals substantial disparities at low kinetic energies, but, as will be 
shown subsequently, bremsstrahlung effects on the general transport process have a significant 
impact only for kinetic energies greater than several MeV. 


2.1.3 Multiple elastic scattering 

Strict conservation principles require that electron-atom interactions generally result in energy 
exchange accompanied by re-direction with respect to spatial variables. Approximations that 
tend to decouple energy loss processes and directional changes greatly simplify electron transport 
analysis. In the present work, these processes are made practically independent. Energy loss is 
specified by collision and radiative losses without regard for directional change, while projectile 
trajectories are described by elastic scattering interactions. Such drastic assumptions call for 
careful scrutiny with regard to both broad energy spectra and the wide variety of material 
types. Some precedents have been set in earlier works [9, 10] in which such approaches have 
been used. In the present formulation, the elastic (Rutherford) scattering cross sections, or,, are 
implemented as a basis for calculating an electron transport mean free path, 


Atr 


P 

na R ’ 


(17) 


where n = pNo/A is the number of scattering centers per unit volume in the medium. See 
Section 3.1 and Ref. [10] for further discussion on transport mean free path. Notice that At r is 
defined to have units of g/cnr. The present development closely follows that of Ref. [8]. The 
screened Rutherford cross section is given as 





Electron Kinetic Energy, MeV 


Figure 2: Radiative stopping power from current model and NIST calculations [6]. 


do R 

dO 


= 2t r 


Ze 2 

( T + m e c 2 )/3 2 ( 1 — cos 6 + 1/e) 


n 2 


where 9 is the scattering angle and e is a parameterized screening function, 

e « 1413[(1 + T/m e c 2 ) 2 - 1 ](Za)“ 2/3 . 

The total cross section is found by integration over all scattering angles. 


or = / — j^(l — cos9)dcos9. 


dO 


-l 


(18) 


(19) 


The integral may be readily evaluated with the substitutions x = cos 9 and K = 1 + 1/e, to give 


-l 


1 — x 
(K - x) 2 


dx = In 


I< + 1 
I< - 1 


1 + 


K - 1 
K + 1 


- ln(2e+1) -fTr 

« ln(2e) — 1 for e > 1. 
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The reciprocal of the transport mean free path, na^/p, may be expressed in terms of kinetic 
energy as 


No N 0 

-OR = 


2irZ z e 


2„4 


[ln(2e) - 1] . 


(20) 


na r _ 

p A R A (T + m e c 2 ) 2 /? 4 

Using expressions from special relativity, terms in the denominator of Eq. (20) may be cast in 
terms of electron kinetic energy. 


/ ? 2 

(T + m e c 2 ) 2 /? 4 


T/m e c 2 , 
T 2 + 2 T 


(T+l) 2 ’ 


(r + l) 4 — 2(r + l) 2 + 1 
(r + 1) 2 

(?n e c 2 ) 4 


(' m e c 2 ) 2 
T 2 + 2 Tm e c 2 — (m e c 2 ) 2 + 


(T + m e c 2 ) 2 ’ 


Now, the electron transport mean free path becomes 


Atr = P/nCTR 


A 


T 2 + 2 Tm e c 2 


2irN 0 e 4 Z 2 


(m e c 2 ) 2 + (m e c 2 ) 4 / (T + m e c 2 ) 2 
ln(2e) - 1 


( 21 ) 


Note that the behavior of the transport mean free path increases monotonically with kinetic 
energy, is nearly logarithmically linear with the log of the kinetic energy, and has generally 
higher values for lower atomic weight elements. 

Another parameter related to elastic multiple scattering is the mass scattering power [11]. 
This quantity is analogous to stopping power, but refers to solid angle scattering rather than 
energy loss. The mass scattering power is defined in Ref. [11] as 


where 


P, = 7T 


2r e Z 

(3 2 (t + 1) 


N 0 

Z 


Y(6r 


( 22 ) 


Q _ Vr 2 + 2t 
r ~ aZ i/3 

with 


Y(e r )^Hi + e 2 )-i + —^. 

The mass scattering power is used in the transport algorithm to relate electron traversal 
on a specific trajectory to the axial penetration distance. The above formulation predicts a 
nearly linear monotonically decreasing variation of the logarithm of scattering power with the 
logarithm of energy. 
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2.2 Photons 


Bremsstrahlung photons produced by electron-atom interactions are attenuated and absorbed in 
the transport medium. In the present work, only the three most important processes involved in 
photon energy degradation and absorption are considered: (1) photoelectric absorption, (2) in- 
coherent scattering, and (3) positron-electron pair production. These processes are dominant for 
different energy ranges, but when taken collectively they closely approximate the total attenua- 
tion for the energy range considered here (0.001 - 1000 MeV). The contribution due to coherent 
(or Rayleigh) scattering has not been included in this analysis since it does not significantly add 
to the total cross section over the energy range of interest (0.001 - 1000 MeV). 

2.2.1 Photoelectric Absorption 

Photoelectric absorption involves interactions of a photon with a bound electron in which the 
photon energy is absorbed with a corresponding kinetic energy increase for the electron. The 
detailed evaluation of the photoelectric cross sections requires the use of complicated atomic 
models in sophisticated solutions of the Schrodinger equation for all energies and elements under 
consideration. For the present formulation, liberal use is made of parametric fitting formulas 
devised from the extensive calculations and tabulations provided by NIST [6]. First, the photon 
energy values for the absorption discontinuities of the K and L shells have been approximated 
as 

Ek ~ 5.435 x 10 _3 Z 2 ' 2038 keV (1 < Z < 92), (23) 

E l « 1.5754 x io~ 4 Z 2 - 5984 keV (30 < Z < 92). 

Note that energies in this section are in keV and L-edge energies are only considered for Z > 30. 
In addition, only the discontinuities due to the K and L shells have been considered. The 
discontinuities (or edges) due to the electronic shells occur at values of the photon energy 
corresponding to the binding energies of electrons involved in the interaction. When the incoming 
photon has an energy greater than the binding energy of a given atomic shell or subshell, a 
discontinuity may occur. 

At the edge discontinuities, the NIST data has also been used to develop formulas for the 
upper (u) and lower (/) values of the photoelectric cross section: 


O'Ku 

= 7165. 5EJ/' 5474 

cm 2 /g, 

(24) 

&KI 

= 593.59-Ej”- 1 ' 327 

cm 2 /g, 


& Lu 

= 9620. 6E^ L5113 

cm 2 /g, 


°Ll 

= 1633.5E^ L2124 

cm 2 /g. 



Examination of typical photoelectric cross sections reveals that the logarithmic slope, 

_ d(ln cr PE ) 
m ~ d { In E v ) ’ 
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is very nearly constant for a given atomic system between absorption edges. Again, using NIST 
cross section data, parametric fits have been used to approximate the logarithmic slopes. 


For E u > Ek: 

/ 6.9068 x W~ 3 Z - 3.014 Z < 10 
mK ~ l — 3.3147Z -0 ' 08156 Z > 11 ' 

For .1 keV < E v < Ek- 


m L 


In 


( a Kl 
y &Lu 



For all Z and Ek < E u < 10000 keV: 

(TpE = (?Ku 



(25) 


(26) 


For Z > 29 and El < E„ < Ek- 

(E v 

(JpE = <JLu I 

An extension to energies less than the lowest edge energy considered is extrapolated with a 
power law function approaching an asymptote based on an extrapolated maximum absorption 
cross section. For Z < 29 and for energies below the A'-edge, an asymptotic cross section is 
calculated using the average ionization potential, 



07 = ®Kl 



The extrapolation formula is then 


(27) 


where the constants are given by 


CPE 


07 

1 + A P E$’ 


(28) 


P = 


rriK 


A P = 


( vi/vki ) - 1 


E p 


(29) 


(VKI/Vl) - 1 

These constants are found by matching the cross section and the slope at the K-e dge. For 
Z > 29, the low energy extrapolation below the L-edge is found by using a L subscript instead 
of K in Eqs. (26)- (29). 

The photoelectric cross sections, as calculated with the present highly parametric formula- 
tion, are shown in Fig. 3 for the three selected materials along with their NIST counterparts. 
Agreement is considered to be fair to good. In circumstances for which low energy photon ex- 
tinction is deemed important, it would probably be best to introduce a tabular database for 
interpolation. 
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Figure 3: Photoelectric absorption cross sections from the current model and NIST calculations [6]. 

2.2.2 Inelastic Photon Interactions 

As the photoelectric absorption coefficient decreases rapidly with increasing energy, the next 
process to dominate is when the photon imparts a part of its energy to an electron resulting 
in scattering of a lower energy photon. The process is referred to as inelastic (or incoherent) 
scattering. The cross section formula for release of an electron of kinetic energy T e by a photon 
of initial energy E v has been taken from Ref. [3] as the free electron Klein-Nishina process. The 
energy constraints are 


0 <T e < T max ; 

_ 2 qE u 
max 1 + 2 q 
E u 

Q — 2 ’ 

where T e is the kinetic energy of the emitted electron and E u is the incident photon energy. The 
cross section formulation is 
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d0KN 

dT e 


NoZ nr 2 } 2 + 


R' = 


A qE u 

NqZ irrl 
A qE u 

T P 


E u - T e 


(E v - T e f 


1 E u — T e 2 (E u — T e ) 

q 2 E v qfT" 


(30) 


B! 


•>\ 2 


R'T P 


2+[ JJ + 


In order to render the cross section applicable to bound atomic electrons, elaborate quantum 
theory calculations are required. Such considerations are only appreciable at lower photon 
energies. For this work, a parametric multiplying function has been developed to approximate 
this effect. 


In Q(Z,E V ) 


0.471 ln(Z) + 1.5184 

K 


[1.003-0. 13317 ln(Z)] 


The final inelastic scattering cross section then becomes 


(31) 


d(Jin 

dT e 


*(Z,E V ) 


dcTKN 

~dT, 


(32) 


For each E u , the total inelastic cross section is found by integration over the appropriate emitted 
electron energy. 


Tmax 

*»= / (33) 

0 

The total inelastic cross sections calculated with the present model are plotted as a function 
of incident photon energy in Fig. 4 along with the corresponding NIST values. The model works 
better for light materials compared to heavy materials and begins to fail at high energies for all 
materials. In this energy range, the inelastic cross section does not contribute appreciably to 
the total photon attenuation, as will be subsequently shown. 

2.2.3 Pair Production 

Electron-positron pair production may occur when a photon of sufficient energy interacts with a 
strong localized Coulomb held. The photon initial energy is transformed into the combined rest 
mass (2 m e c 2 ) of the newly created particles. Any remaining energy appears as kinetic energy 
of the new particles. A comprehensive description of the process may be found in Ref. [12], 
where a variety of cross section quantities are derived for several aspects of the phenomenon. 
The formulation from Ref. [11], chosen for use in the present work, applies to the cross section 
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Photon Energy, MeV 


Figure 4: Inelastic photon cross sections from the current model and NIST calculations [6]. 


for production of a positron of total energy E + by a photon of initial energy E u . Only positron 
production associated with pair production is considered in this work. The pair production 
differential cross section may be written in simple terms as 


where 


do 'pp 

~dE + 


N 0 ar 2 e Z(Z + Q 

A El 


(Q i + Q 2), 


(34) 


Qi = [E\ + (E u - E+fW^x - F c ], 

Q2 = \e + {E v - T + )[02 - F c \. 

The terms £, 4> 2 ( 8 ') and F c are the same as for the bremsstrahlung cross section, Eqs. 

(11)-(13), with the exception that the value of S' is determined as 


136 m e c 2 E v 
Z 1 /‘ i E + (E u — E + ) ' 


Energy conservation constrains the positron energy as 

TYleC 4 E/_(_ T T max -)- Tl~leC , 


T 

^ max 


= E u — 2 m e c 2 . 


(35) 
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The above formulas for pair production become inaccurate as the pair production threshold 
is approached. Consequently, a power law function is fit at the first energy grid value greater 
than 4 MeV (= £ 4 ). Now, the low energy pair production cross section is given as 

app = A pp (E u - 1.03) 6 for 1.03 <E V < £ 4 , (36) 

with 

da pp £4 — 1.03 

dE+ e+=Ea cr pp (£4) 

<Tpp(£ 4 ) 

pp (£4 — 1.03) b ' 

A source term for positron production may be derived in a manner similar to that described 
previously for production of bremsstrahlung photons. See Eq. (15). For a photon differential 
flux spectrum (p v {E v ,x), positron production may occur between energies £ pP) min = 2 m e c 2 and 
£pp,max = £^,max — 2 m e c 2 , where £j, imax is the highest photon energy value considered. The 
expression for the source term becomes 

£pp,max 

C (£+,*)= J Vv{Kix)jj^dE'„. (37) 

^"pjp . rn i r i 

It should be noted that in parallel with the positron production process, a similar cross 
section expression applies to the creation of the partner electron, along with the corresponding 
source term for secondary electrons. Just as for the inelastic process, the pair production cross 
section may be found by integration over the allowed positron energy ranges 

E+,max 

»pp =j ^ dE' + , (38) 

m e c 2 

where £+, m ax = E v - m e c 2 . 

Results of pair production cross sections, as calculated from the present model, are shown in 
Fig. 5 along with the appropriate NIST comparison values. Agreement is very good except for 
energies near the threshold region where the power law extrapolation formula has been invoked. 

2.3 Cross Section Implementation 

The individual cross sections described previously are used in combination to provide three crit- 
ical parameters essential for description of the electron-photon transport process. The ultimate 
slowing and stopping of electrons is governed by the total stopping power obtained from the 
sum of the collision and radiative stopping powers. The collision process dominates at low en- 
ergies (less than approximately 1 MeV), while the radiative process assumes the dominant role 
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Figure 5: Pair production cross sections from current model and NIST calculations [6]. 


at high energies (greater than approximately 1 MeV). A transitional minimum occurs in the 
neighborhood of 1 MeV as is shown in Fig. 6, where the total stopping power is plotted along 
with the comparison values taken from Ref. [6]. Note that the lower energy inaccuracies of the 
model radiative stopping power have insignificant impact on the total stopping power. 

The general attenuation and extinction of photons in a medium is found to be closely re- 
lated to the sum of the photoelectric, inelastic, and pair production cross sections. Coherent 
scattering processes have not been considered because of their relative unimportance in space 
radiation effects. Fig. 7 shows the model total attenuation coefficient (taken to be the sum of the 
photoelectric, inelastic and pair production cross sections) and the selected comparison data, 
where the model results are again deemed satisfactory for an adequate description of the photon 
transport. 

The final parameter of importance to the electron-photon transport process is that of the 
photon energy deposition coefficient which is used to calculate effective photon dose in Eq. (64). 
This relates to those fractions of the photon attenuation processes (fpE, fin, / pp ) that produce 
secondary electrons. For photoelectric absorption, it is assumed that the secondary electron is 
emitted with the same energy as the incident photon (i.e. the binding energy of the released 
electron is neglected) and the value of /pe is unity. In the case of inelastic scattering [13], 


/in = 1 - 


E, 


v, in 


E„ 


where E u ; n is the energy of the photon produced in the inelastic process initiated by a photon 
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Electron Kinetic Energy, MeV 


Figure 6: Model total stopping powers compared with NIST calculations [6]. 


of energy E v . For the pair production process, / pp = 1 — 2 m e c 2 /E u . 

An additional reduction in the photon energy deposition coefficient is required to discount 
those secondary electrons that go on to produce additional bremsstrahlung. In the present 
model, this factor is represented by (1 — Srad/Srot)) w here Sr 0 t = Srad + <Scoi • This factor is 
only applied to the inelastic and pair production processes, so that the model expression for the 
total photon energy deposition coefficient becomes [13] 

Tnax.in Tmax,pp 

M.» = VPE+J (1 - + / (1 - f^)/pp^r', (39) 

0 0 

where the upper limits on the integrations apply to the maximum allowed emitted electron 
energies for the respective processes. The energy deposition coefficient is described in detail in 
the X-Ray Data section of the NIST website [6]. A comparison of model results with the more 
elaborate NIST calculations is shown in Fig. 8. 

3 Transport Formulation 

In the present work, the general transport process is focused principally on description of pen- 
etration of a primary electron field along with secondary bremsstrahlung photons generated by 
electron-atom interactions. For electrons, an essentially one-dimensional formulation is devel- 
oped with some reference given to the effects of radially symmetric scattering. While the electron 
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Photon Energy, MeV 

Figure 7: Model photon total attenuation coefficient compared with NIST values [6]. 


propagation is determined by the stopping powers, the photons are assumed to be transported 
along the direction of travel of the electrons. Photon intensity is governed by the radiative 
transfer equation [14], which utilizes the calculated attenuation coefficients. 


3.1 Electrons 

A quantity of fundamental importance in describing transport of electrons in matter is that of 
the maximum distance of travel as determined from the energy loss stopping powers calculated 
from the modified Bethe-Bloch and bremsstrahlung formulas in Eqs. (1) and (16). This quantity 
is referred to as the Continuous Slowing Down Approximation (CSDA) range R(E). For a given 
energy, E, the CSDA range is determined from [2] 



/ 


dEf 


*Scol + Srad 


(40) 


The inverse of the R(E) function, i?” 1 (E) = E(R), may be used to construct the variation 
of dE/dx as a function of distance traveled by an electron of initial energy E$. Such a function, 
S(Eq,x), is analogous to the usual Bragg curve, but neglects straggling and measures energy 
loss as a function of distance traveled along the particle’s path rather than perpendicular depth 
in a material. Straggling is accounted for after the incorporation of the statistics associated with 
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Figure 8: Photon energy deposition coefficients from model compared with NIST values [6]. 


scattering. The residual energy of an electron, at position s, is given by 


W(s) 


R 

j S(Eo,x)dx 
o 


J S(Eo,x)dx = Eq — 
o 


S 

J S(Eo,x)dx, 
o 


s < R. 


(41) 


For a beam of mono-energetic electrons, it is assumed that the maximum distance traversed 
is the CSDA range and straggling is neglected. However, the average distance of penetration 
along the beam axis direction is generally less than the CSDA range due to multiple scattering 
effects. The average deflection for a unit path length can be defined by the inverse mean free 
path (A,^ 1 ) [9] which yields 


and 


(cos 0) 


d{z) 

ds 


= exp 


S 


0 


ds 

Atr {s') 


(z(s)) 




ds 


(42) 


(43) 


In general, (x) is defined to be the arithmetic average of the variable x. 


For electrons, ( z(R )) 
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is interpreted as an “effective range” or the axial distance at which 50% of the electrons have 
stopped. In other words, the z-value for which the transmission is 0.5. 

A further critical assumption is made by specifying a Gaussian distribution about ( z(R )), 
where ( z(R )) is the distance of peak electron population. In addition, the condition that practi- 
cally all electrons are stopped at z = R(E ) is assumed. This Gaussian represents the variation 
in electron path length due to multiple scattering effects. By invoking the formula for half-width 
of the Gaussian distribution, a value for deviation, H, is found as [15] 

^_ R-(z(R)) 

2v/2 In 2 

The resultant Gaussian function may be interpreted as the probability that an electron pene- 
trates an axial distance z within increment dz, 


P{z) 



liiM 

2K 2 


(45) 


Integration of the Gaussian function results in an error function prescribing the number of 
particles having stopped over distance z. The complementary error function [16] results in the 
corresponding transmission function, 


r](z) = 0.5 < 1 — erf 


( Z ( R )) 

K \/2 

A special algorithm for the standard error function has been utilized in the present code, 
been extracted directly from Winitzki [17]. 


(46) 
It has 


erf(x) 


1 — exp 



a x 2 + £ 

ax 2 + 1 


1 1/2 


where 


8 7T — 3 
37T 4 — 7T 

For the negative argument, the symmetry relation erf(-.x) = -erf(x) is invoked. Sample trans- 
mission functions have been calculated for mono-energetic electrons incident on aluminum and 
are shown in Fig. 9. The results indicate the marked increase of (z) as a fraction of R with 
increasing initial energy. 

In addition to specifying the transmission probability of electrons as a function of axial 
penetration distance z, it is necessary to evaluate the variation of primary electron energy as 
a function of 2 . A formula has been developed [10] in terms of the moments of the mass 
scattering power that expresses, in cylindrical coordinates ( z,r ), the distance covered for an 
average electron trajectory having axial distance z and deflection radius r. The moments are 
calculated as 
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z/R 


Figure 9: Electron transmission functions in aluminum calculated for three initial energies. 


a i = J P s {u)(z — u) l du, 
0 


(47) 


with P s given by Eq. (22). The moments are then used to derive an equation, details of the 
derivation are given in Ref. [10], for distance along an average trajectory passing through ( z,r ), 


s(z , r ) 


2: 


0 


1 + ao(z ) + 



a i( z ') 

a2(z') 


1 1/2 


dz . 


(48) 


The quantity \fa?[z) has been identified [10] as the half- width at half- height of the distribution 
of deflection radius at axial distance z. This observation leads to the inference that a “most 
typical” trajectory intersects the point (z, ydx^). The third term in Eq. (48) then vanishes and 
an “average mean path” becomes, 


(s(z)) = J \Jl + a 0 (z)dz. (49) 

o 

Several of these “typical mean trajectories” have been calculated for a range of incident 
kinetic energies for the three materials chosen for illustration. The plots of Figs. 10, 11, 12 
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Figure 10: Typical mean trajectories for electrons in H 2 O. 
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Figure 11: Typical mean trajectories for electrons in aluminum. 
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Figure 12: Typical mean trajectories for electrons in tantalum. 


display the quantity y/a^ versus axial distance z. The curvilinear distance (s(z)) along each 
path is then represented by Eq. (49). 

Reference to Eqs. (43) and (49) shows that the relationship between the mean axial distance 
and average mean path is such that (z) < z < (s). The average residual energy (W) may be 
associated with an average distance along a path length ( s ) by again interpolating on a table of 
CSDA range-energy relations. However, the interpolated value (W) reaches zero before z = R 
is attained since (s) is a value that includes average deflection and z = R is the axial distance 
traveled in the continuous slowing down approximation assuming no angular deflection. The 
axial penetration, however, may extend beyond the distance at which (W) vanishes due to energy 
straggling. A provision for energy straggling is supplied by activating a power law function of 
the form 


(W) = A W (R- z) k . (50) 

This formula is invoked somewhat arbitrarily for z - values for which (W) < 0.5 Res DA • When 
this condition pertains, the constants A w and k are evaluated by point-slope matching at the 
z-value where (W) = 0.5 R^csda- Fig. 13 shows sample normalized residual energy curves 
calculated for electrons in aluminum at the initial energies of 5 and 50 MeV. The 5 MeV results 
show considerable divergence of (W) from the CSDA loss curve, whereas (W) for 50 MeV adheres 
much more closely to the corresponding CSDA function. The transition to the energy straggling 
extension is apparent for both cases. 

The above transport parameters may be used to devise an expression for the differential flux 
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Figure 13: Residual energy functions for electrons in aluminum for 5 and 50 MeV. 


spectrum, t p(z,T ), of electrons at distance z, given the initial spectrum at z = 0, ip(0,To). In 
going from the boundary to distance z, particle number conservation requires that (p( 0, To)<5To = 
p>(z, T)5T, where the energy scaling corresponds to that prescribed by the CSDA range-energy 
relation for the material. In general, SY is defined as a small, finite difference in the value of 
Y . In addition, the transmission function may be used to account for particles not arriving at 
z due to interaction not accounted for in the stopping power. 

?(*.r ) = MWvMmW (si) 

Division of the kinetic energy increments by distance increments leads to the ratio of stopping 
powers in Eq. (51). Thus, for a slab of thickness w, the differential spectrum of primary electrons 
may be found at any axial distance z < w. If primary electrons were the only concern, the 
calculation could end with Eq. (51). However, bremsstrahlung photons are generated during the 
transport process. In order to take this process into account, a spatial grid is established within 
the slab layer. In the present code, a slab layer is assigned a spatial grid, normally of 20 to 30 
points, with spacing increasing monotonically according to the following formula, 

d ‘ = (wz t ) ”• 

Here, di is the depth in the material slab of the zth point and N is the total number of points 
in the spatial grid of the slab. The photon source term may then be calculated at each spatial 
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grid point as 


T{di) 

C (E„,di)= J ipeiT'j^dT for 0 <E v <T(di). 

E v 

Finally, the primary electron spectra provided by the present code are calculated at each 
spatial grid point as 

fJ ^ _ (p(0,Eo)rj(di,Eo)S{Eo) 

<Pe,i{dt,T) — S(T) 

Note that the transmission function r/ has been defined for each initial energy grid value, with 
the kinetic energy T(di) being the corresponding residual energy at d t . Currently, the code 
accounts for secondary electron production in the dose calculations by assuming the electron is 
deposited locally at the point of production. 

Implementation of Eq. (51) for the illustrative LEO spectrum [18] normally incident on 
aluminum provides the results shown in Fig. 14 (the environment spectrum is labeled Z = 
0.0). The stopping of low energy electrons and the transfer of high energy electrons to lower 
energies is evident in the spectral functions at increasing penetration depths. Similar results for 
the Europa spectrum in aluminum are shown in Fig. 15 using an environment from the Galileo 
Interim Radiation Electron (GIRE) Model [19]. The much greater penetration and more gradual 
attenuation is clearly evident for this very high energy environment. 

3.2 Photons 

In the previous section on cross section evaluation, it was shown that a spectral source term 
for bremsstrahlung photons may be calculated, Eq. (15). When the differential electron flux 
spectra have been specified as a function of axial distance z, a corresponding photon flux may 
be determined using a solution of the radiative transfer equation for an emitting-absorbing 
medium [14]. In addition to the source function, the total attenuation coefficient is required. 
The total attenuation coefficient is the sum of the relavant cross sections, 

C’Tot = ^PE T °"in T <7pP' (52) 

The photon flux spectrum at each spatial grid point di then involves integration of the source 
term multiplied by the attenuation factor. 

di 

<Pv,i(di, E v ) = J ({x, E v ) exp [-OTot {<k - a;)] dx. (53) 

o 

Photon spectra for the Europa spectrum are illustrated in Fig. 16 and correspond to the 
electron spectra of the preceding figure. The persistence of higher energy photons at large 
aluminum depths is indicative of their highly penetrating nature. 
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Flux, # electrons/(cm 2 MeV day) Flux, # electrons/(cm 2 MeV day) 



Kinetic Energy, MeV 

Figure 14: LEO electron spectra at several aluminum thickness values. 



Figure 15: Europa electron spectra at several aluminum thickness values. 
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3.3 Positrons 


Ordinarily, exposure effects due to space environment electrons may be adequately evaluated 
from specified electron and photon flux values at a given target location. The pair production 
contribution to the energy absorption coefficient should provide a fair estimate of exposure due to 
positrons. However, if details of specific positron interactions are of interest, it becomes necessary 
to include explicit transport of the positrons and their associated annihilation photons. In the 
present code, source terms and flux for positrons are evaluated, along with source terms for the 
annihilation photons, in an uncoupled manner with a view toward allowing the user to assess 
the relative importance of the presence of positrons. The positron source term, C(E+, x), has 
already been given as Eq. (37). These antimatter particles have a “sink” term by annihilation 
with the constituent bound electrons of the medium. This term may be represented as [8] 


N e cr ann (r ) 


N e nrl 
t + 2 


t 2 + 6r + 6 
t(t + 2) 



(54) 


where N e is the number of electrons per unit mass of material and cr ann (r) is the total annihilation 
cross section, with 


h\ = In 


T + 1 + \/t{t + 2) 


h 2 


T + 4 

\J t(t + 2) ’ 


and t representing the positron kinetic energy in rest mass units (0.511 MeV). 

Positrons initially at spatial location Xj t arriving at xi + 5x with energy Ej began traverse of 
this increment at a higher energy, Eq = E{R + dx), where R is the residual range at x t + 5x. 
The CSDA process over the interval 5x with no loss of positrons requires that 


ip(Eo,Xi)S(E 0 ) = p(Ej,Xi+i)S(Ej). 


(55) 


When annihilation takes place, the fraction of positrons from the incident spectrum lost in 
distance Sx is 


At(5x) = exp 



The positrons transmitted through interval 5x is then 


c p^(Ej,Xi + Sx) = psA T (5x)(p [E ( Rj + Sx ) , Xi ] , 


(56) 


(57) 


with 


S [E ( Rj + Ax)] 

S(E~) 
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Additional contributions to the positron flux at Xi + \ = Xi + 5x arise from the distributed 
sources in the interval that must include an additional annihilation factor, Ag. 


Sx 

(Ej,x + Sx) = j C[E(R + x ),Xi + x ]Ag(x )dx , (58) 

o 

with 



/ r 

X 


' 

Ag(x ) = exp < 

/ N e <7 ann 

E (^Rj + x ^ 

dx > 


J 

0 

V 




Here, Ag is the fraction of positrons lost from the positrons created in the distance increment. 

The total positron flux, ip e + = tp® + <p^ s \ is calculated in the code by simple numerical 
integration procedures applied to each of the established spatial grid intervals. Results for the 
Europa spectrum for selected thicknesses in aluminum are shown in Fig. 17 and may be compared 
with the corresponding photon flux of Fig. 16. It is seen that the positron flux is a relatively 
small fraction of the governing photon field and that the uncoupled calculation is justified. 

When the positrons are annihilated, photons are generated in accordance with the differential 
cross section [8]. 


where 


^ann 

dE u 


* r l 

t(t + 2) 


[?(«) +q(r + 2 




q(x) 


1 

x 


+ 2 + 2 


T + 1 
r + 2 



(60) 


with r and k being the positron kinetic energy and the photon energy, respectively, in units of 
electron rest mass. The annihilation photons are restricted by kinematics to an energy range of 


m e 

1 + a' 


< E 


v, ann 


< 


m e 

1 -a r 


(61) 


where 


a 


I 


r + 2 


Finally, a source term for the annihilation photons per unit mass of material may be expressed 
as 
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Flux, # positrons/(cm 2 MeV day) S Flux, # photons/(cm 2 MeV day) 



0.01 0.1 1 10 100 

Photon Energy, MeV 


16: Photon spectra from Europa electrons [19] at several aluminum thickness values. 



Kinetic Energy, MeV 


Figure 17: Positron flux at selected thicknesses in aluminum for Europa spectrum [19]. 
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Tmax 

C,ann (E v ,x) = f p(E + , x) ^dT + . (62) 

T • 

I min 

Production of annihilation photons occurs with highest probability near the electron/positron 
rest mass energy. For the Europa case of Fig. 17, the annihilation photon source terms calculated 
according to Eq. (62) are shown in Fig. 18 and exhibit the expected peak at 0.511 MeV. Since 
this contribution is small relative to the other processes, these photons are not transported in 
the present formulation. 

4 Dosimetric Calculations and Comparisons 

The conventional radiation dose for ionizing radiations in matter is defined as the energy im- 
parted to a mass element of the material by the particles traversing the elemental mass. In 
general, only energetic charged particles are responsible for the energy transfer. The “dose” 
attributed to uncharged particles (e.g. photons, neutrons) results from the charged particles 
generated within the medium by the neutrals. 


4.1 Electron and Photon Dose Evaluation 

Energy deposition of electrons is calculated by multiplying the local flux (differential in energy) 
by the total stopping power and integrating over energy. The dose at a given location is given 

by 


oo 

D e (x) = J ip e (x,E)S(E)dE. (63) 

o 

In Eq. (63), E represents the electron energy at position x. The local differential flux, <p e , has 
units of electrons/ (cm 2 MeV). The stopping power, S, has units of MeV cm 2 /g, while the units 
of D e are MeV /g which may be converted to cGy or rads upon multiplication by the conversion 
factor 1.602 x ICG 8 . For the photons, absorbed dose is characterized by the energy deposition 
coefficient, /i en , which is described in Eq. (39) and is expressed herein as a mass absorption 
coefficient with units cm 2 /g. The dose attributed to photons is then 

OO 

D u (x) = J Hen E u <p v (E v , x)dE u . (64) 

o 

The use of Hen to calculate effective photon dose (as in the present version of the code) is a 
simplistic approximation method that has been widely used and considered to be an adequate 
representation of this dose contribution [3,6]. This approximation assumes that all the energy 
deposited by photons is deposited by charged particles at the point of their production. A more 
direct evaluation would be to use the appropriate cross sections for generation of secondary 
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Figure 18: Annihilation photon source for Europa spectrum [19] in aluminum. 


electrons (ope, Om, and cr pp ) to compute respective source terms in the manner expressed in 
Eq. (37). The corresponding flux terms could then be determined using procedures analogous 
to that described for the positron transport in Eqs. (54)- (58). Such modifications may be con- 
sidered for future upgrades if the additional complexity appears to be warranted. 


4.2 Sample Calculations and Comparisons 

Calculations performed with the present code were selected to examine the behavior of trapped 
electron spectra interacting with light and heavy materials. The chosen electron environment 
spectra (boundary conditions) are shown in Figs. 14 and 15 and represent, respectively, the 
environment in LEO (400 km, 51° inclination) [18] and the environment in the vicinity of the 
orbit of the Jovian moon Europa [19]. The initial environment is shown as the curves labeled 
Z=0.0 in Figs. 14 and 15, along with the flux spectra at various depths in aluminum. After 
transport through the shield medium, the doses have been evaluated in silicon to simulate 
exposure of a solid state device. 

The Integrated TIGER Series (ITS) Monte Carlo code [20] was used to calculate dose for 
the LEO spectrum electrons on aluminum at normal incidence for several depths. The identical 
scenario was used with the LaRC deterministic code and the results are compared in Fig. 19. The 
LaRC code gives generally lower values for the dose versus depth functions than the Monte Carlo 
results, but the functional behavior is very similar. Greater differences are to be expected for 
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the low energy spectra of LEO because scattering processes are more prominent in the transport 
process and are treated differently in the two calculations. 

A second exercise using the LEO spectrum was performed using the SPENVIS space envi- 
ronment website that incorporates the application code SHIELDOSE-2 that is based on results 
from ITS, parameterized for several simple geometries [18]. The geometry selected is that of 
the semi- infinite slab adjacent to an isotropic field of energetic electrons having the same LEO 
energy spectrum previously used. For the LaRC comparison calculation, a dose versus depth 
function was generated for normal incidence on aluminum. A solid angle integration was then 
performed over a directional hemisphere about a target point on the shielded side of the slab. 
For a slab of thickness w, the calculation may be expressed as 

2n 7 r 1 

D( - w) = hJ J D (^) s ™ md * = lJ D (^d) dcosl> - (65) 

0 0 0 

The comparison of the LaRC results with those of SPENVIS are shown in Fig. 20 and again 
indicate that the LaRC results for total dose are generally lower than those based on the Monte 
Carlo calculation. This adds credibility to the previous normal incidence calculations of Fig. 19. 

For higher energy electron spectra, sample application comparisons have been made for 
a Jovian electron environment generated by the NASA-JPL GIRE model [19]. The relevant 
spectrum is shown in Fig. 15 as the curve labeled Z = 0.0 and represents electron flux in 
Jupiter’s equatorial plane at the average orbital distance of Europa. The normal incidence 
on a semi-infinite slab scenario was used for the materials aluminum, tantalum, and copper- 
tungsten (50%-50%) alloy. The previous LEO spectrum exhibited very few electrons above 5 
MeV, whereas the Europa electrons may have substantial population up to and beyond 100 
MeV. 

The dose versus depth curves for the Europa spectrum are given in Figs. 21, 22, and 23 for the 
specified materials. In general, the LaRC deterministic calculations show improved agreement 
with the corresponding Monte Carlo results compared to the lower energy LEO environment. 
The range of scaled thickness for the Europa cases is from 0.5 to 20 g/cm 2 . It is seen that 
electron stopping is practically complete at the end of this thickness range, after which the 
bremsstrahlung contribution dominates. 

5 Concluding Remarks 

The electron/photon transport code described here has been developed with a view toward 
simplicity and speed in analysis of exposure from space environment electrons. In its present 
form, the code may be implemented to great advantage in shield material trade studies, numerical 
and statistical experiments, uncertainty analyses, etc., despite some observed differences with 
Monte Carlo comparisons. Numerous comparative calculations, in addition to those reported 
here, have been performed with similar degrees of agreement as those presented in the text. 
In particular, calculations for the trapped Jovian environment compare more favorably with 
corresponding Monte Carlo results than do the comparisons for the much lower energy-range 
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Figure 19: Comparison of dose rate results in aluminum for the LaRC model with the ITS [20] model 
from a LEO spectrum [18]. 



Figure 20: Comparison of dose rate results in aluminum for the LaRC model with the SPENVIS [18] 
model from a LEO spectrum [18]. 
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Figure 21: Comparison of dose rate results in aluminum for the LaRC model with the ITS [20] for 
Europa spectrum [19]. 



Figure 22: Comparison of dose rate results in tantalum for the LaRC model with the ITS [20] for Europa 
spectrum [19]. 


35 



Material Amount, g/cm 2 


Figure 23: Comparison of dose rate results in W-Cu alloy for the LaRC model with the ITS [20] for 
Europa spectrum [19]. 


spectra of LEO. This fact suggests that the deterministic formulas for low-energy scattering 
should be improved in future versions. 

The present code may be used in a variety of space mission applications, but may also be 
considered as developmental for several reasons. The various cross section comparisons with 
NIST calculations show several instances where the LaRC code cross sections may be improved 
by introducing selective correction terms, even though the cross section deviations noted have 
little impact on final exposure results. With regard to transport, the present code is amenable 
to immediate extension to a two-dimensional axi-symmetric representation. In addition, more 
explicit details pertaining to very high energy processes involving positrons and their annihilation 
photons may be obtained by a formal coupling of these species to the general electron/photon 
transport. Although no significant errors are indicated as a result of the decoupling of positron 
and annihilation photon transport for the Jovian spectral environments, use of the present code 
for high energy beam simulation or cosmic ray shower analysis would most likely require full 
coupling of these processes. It is therefore natural to anticipate code upgrades and extensions 
in the near future. 
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Appendix A Input File of Atomic Data 


Table 2: Input File of Atomic Data used in the Electron Transport code. Values appearing in 
italics for the material density, p, indicate the value corresponds to the condensed phase of the 
material. 


Sym 

Z 

A 

P (g/cm 3 ) 

I (eV) 


Sym 

Z 

A 

p (g/cm 3 ) 

I (eV) 

H 

1 

1.008 

0.071 

19.2 


Ag 

47 

107.868 

10.5 

470 

He 

2 

4.003 

0.122 

41.8 


Cd 

48 

112.41 

8.65 

469 

Li 

3 

6.941 

0.534 

40 


In 

49 

114.82 

7.31 

488 

Be 

4 

9.012 

1.848 

63.7 


Sn 

50 

118.71 

5.75 

488 

B 

5 

10.81 

2.37 

76 


Sb 

51 

121.75 

6.691 

487 

C 

6 

12.011 

2 

78 


Te 

52 

127.6 

6.24 

485 

N 

7 

14.007 

0.808 

82 


I 

53 

126.905 

4.93 

491 

O 

8 

15.999 

1.14 

95 


Xe 

54 

131.29 

3.52 

482 

F 

9 

18.998 

1.108 

115 


Cs 

55 

132.905 

1.873 

488 

Ne 

10 

20.179 

1.207 

137 


Ba 

56 

137.33 

3.5 

491 

Na 

11 

22.99 

0.971 

149 


La 

57 

138.906 

6.145 

501 

Mg 

12 

24.305 

1.738 

156 


Ce 

58 

140.12 

6.657 

523 

A1 

13 

26.982 

2.699 

166 


Pr 

59 

140.98 

6.7 

535 

Si 

14 

28.086 

2.33 

173 


Nd 

60 

144.24 

6.9 

546 

P 

15 

30.974 

2 

173 


Pm 

61 

145 

7.22 

560 

S 

16 

32.06 

2 

180 


Sm 

62 

150.36 

7.5 

574 

Cl 

17 

35.453 

1.56 

174 


Eu 

63 

151.96 

5.243 

580 

Ar 

18 

39.948 

1 

188 


Gd 

64 

157.25 

7.9 

591 

K 

19 

39.098 

0.862 

190 


Tb 

65 

158.925 

8.229 

614 

Ca 

20 

40.08 

1.55 

191 


Dy 

66 

162.5 

8.55 

628 

Sc 

21 

44.956 

2.989 

216 


Ho 

67 

164.93 

8.795 

650 

Ti 

22 

47.88 

4.54 

233 


Er 

68 

167.26 

9.066 

658 

V 

23 

50.942 

6.11 

245 


Tm 

69 

168.934 

9.321 

674 

Cr 

24 

51.996 

7.19 

257 


Yb 

70 

173.04 

6.9 

684 

Mn 

25 

54.938 

7.3 

272 


Lu 

71 

174.967 

9.84 

694 

Fe 

26 

55.847 

7.874 

286 


Hf 

72 

178.49 

13.31 

705 

Co 

27 

58.932 

8.9 

297 


Ta 

73 

180.948 

16.654 

718 

Ni 

28 

58.69 

8.902 

311 


W 

74 

183.85 

19.3 

727 

Cu 

29 

63.546 

8.96 

322 


Re 

75 

186.207 

21.02 

736 

Zn 

30 

65.39 

7.133 

330 


Os 

76 

190.2 

22.57 

746 

Ga 

31 

69.72 

5.904 

334 


Ir 

77 

192.22 

22.42 

757 

Ge 

32 

72.59 

5.323 

350 


Pt 

78 

195.08 

21.45 

790 

As 

33 

74.922 

1.97 

347 


Au 

79 

196.967 

19.3 

790 

Se 

34 

78.96 

4.79 

348 


Hg 

80 

200.59 

13.546 

800 

Br 

35 

79.904 

3.12 

357 


Tl 

81 

204.383 

11.85 

810 

Kr 

36 

83.8 

1 

352 


Pb 

82 

207.2 

11.35 

823 

Rb 

37 

85.468 

1.532 

363 


Bi 

83 

208.98 

9.747 

823 

Sr 

38 

87.62 

2.54 

366 


Po 

84 

209 

9.32 

830 

Y 

39 

88.906 

4.469 

379 


At 

85 

210 

1 

825 

Zr 

40 

91.224 

6.506 

393 


Rn 

86 

222 

4.4 

794 

Nb 

41 

92.906 

8.57 

417 


Fr 

87 

223 

1 

827 

Mo 

42 

95.94 

10.22 

424 


Ra 

88 

226.025 

5 

826 

Tc 

43 

98 

11.5 

428 


Ac 

89 

227.028 

10.07 

841 

Ru 

44 

101.07 

12.41 

441 


Th 

90 

232.038 

11.72 

847 

Rh 

45 

102.906 

12.41 

449 


Pa 

91 

231.036 

15.37 

878 

Pd 

46 

106.42 

12.02 

470 


U 

92 

238.029 

18.95 

890 
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